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Abstract 

We present time-dependent analytic solutions to the Einstein equa¬ 
tions coupled with a dilaton (scalar) field. The background geometry 
for the solutions is a product of an A-dimensional spherically symmetric 
space and a d-dimensional flat space. We discuss the global properties of 
the spacetime. 


1 Introduction 

In recent years various types of inhomogeneous cosmological solutions to the 
Einstein equations coupled with other fields have been studied extensively [1, 2]. 
Such exact solutions are important for studying the structure and properties of 
horizons, because knowledge of the global structure of the spacetime is essential 
for such a study. 

On the other hand, time dependent solutions to the Einstein equations cou¬ 
pled to matter fields have been investigated as models of gravitational collapse 
[3]. A convenient model may involve a scalar field as a matter field, which is 
coupled to Einstein gravity. The gravitation theory including scalar fields is pre¬ 
ferred not only because of its simplicity, but also because it can be considered 
as a reduced system of supergravity theory [4] or superstring theory [5]. 

Exact solutions to the coupled Einstein massless scalar field equations have 
been obtained for some simple cases, including a static case [6]. Besides a 
homogeneous cosmological solution, only a few time-dependent exact solutions 
describing an inhomogeneous spatial metric are known. The first example of 
such a solution has been given by Roberts [7]. Another type of exact solution 
has recently been given by Husain, Martinez and Niihez [8]. 
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In this paper we will obtain the time-dependent spherically symmetric exact 
solution to the Einstein scalar theory in arbitrary spacetime dimensions. This 
type of solution is a generalization of the one found by Husain, Martinez and 
Nunez [8].^ We will consider the metric in Kaluza-Klein pararmetrization [9], 
which represents topologically a product space of a (1 -b iV)-dimensional space 
and a d-dimensional flat space. The configuration of the scalar field and the 
metric is assumed to be spherical in the N dimensional subspace. 

As well as the massless case, in this paper we consider the case with a 
potential term represented by an exponential function of the scalar field variable. 
Such a potential often arises from effective field theories of string theory or 
supergravity theory. In these theories the scalar field is known as a dilaton field 
[4, 5]. The exponential potential plays a crucial role in the exact multi-centred 
solution [2] as well as in exact solutions for cosmological inflation [10]. 

One motivation for studying the time dependent solutions to the coupled 
Einstein scalar system is that they may be regarded as an analytic model for 
gravitational collapse induced by scalar fields. Recently there has been much 
progress in the numerical study of gravitational collapse. Numerical results 
illustrate that the behaviour of a scalar field configuration may have a certain 
self-similarity and the critical mass of the resulting black hole may take a form 
governed by a certain power law with a universal critical exponent [11]. Several 
authors have made efforts to understand the critical phenomena qualitatively 
by using the exact solutions. Brady [12] and Oshiro et al [13] have discussed the 
critcal exponent, based on he exact solution given by Roberts [7], while Husain 
et al based their discussion on the exact solution derived by themselves [8]. 

Throughout this paper we consider the action, including a dilaton field </>, of 
the form: 


I = 


d^: 


.yEl 

1071 


R- 


D-2^ 


-b boundary terms 


( 1 ) 


where D = \ N d, and the Newton constant is set to unity. Here we assume 
that the dilaton coupling a may take an arbitrary positive value: for effective 
field theories of string theory it takes a = 1. Eor a = 0, A simply becomes a 
cosmological constant. 

In the following section we first consider the A = 0 case, where the dilaton 
field ^ is a free massless field. We will treat the case with A ^ 0 in section 
3. In section 4 we discuss the structure of singularities and apparent horizons 
in the spacetime described by the exact solutions obtained in sections 2 and 3. 
In section 5 we evaluate the mass for the self-gravitating system described by 
exact solutions with the massless scalar field obtained in section 2. 

^It is difficult to generalize Roberts’ solution to the arbitrary dimensional case. 
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2 Solutions for a massless scalar field 


For A = 0 the field equations derived from the action (1) take the following 
form: 


□ ^ = 0 

4 

Rmn = • 


( 2 ) 

(3) 


We wish to find the time-dependent solution to the equations (2) and (3) 
which can be interpreted as a product of an A^-dimensional spherically symmet¬ 
ric space and a d-dimensional flat space.^ We assume that the metric should 
take a block diagonal form in d- and A^-dimensional spaces. Throughout this 
paper we use the following metric; 


ds^ = \-/\{r)dt^+T^{t)dx^]W/^^~^\r) S^(t) 


dr"^ 

_A(r) 




(4) 

where dx^ = dx^, and dn‘jf_^ represents the line element of a unit (N—l)- 
sphere. Here the scale factors S and T are assumed to be functions of t, while 
A and cr are functions of r. 

Now we can transform the equations (2) and (3) into simultaneous differen¬ 
tial equations concerning the unknown functions S, T, A, cr and (p, by taking the 
metric ansatz (4). Furthermore, we require the separation of equations with the 
variables t and r to obtain analytic solutions. To this end we will take another 
ansatz for the scalar field: the time derivative of </> depends only on t, while the 
r-derivative of (p depends only on r. We also note that the equations for the 
scale factors should coincide with those known in the Kaluza-Klein cosmological 
scenario [4, 9]. 

Consequently, the reduced equations we obtain are divided into three groups. 
One group includes the equations for the temporal evolution of S and T: 


1 


gN'J'd 


= 0 


N- 


S 


T 4 

~ ~ 0 - 2 ^’ 


I 




T 

T 


-I- (d - 1) 


T 

T 


ST 


where a dot denotes derivation with respect to t. 


^For static cases several exact solutions are obtained in [14]. 


(5) 

( 6 ) 

(7) 

( 8 ) 
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Another group includes the equations for the spatial configuration of A and 


a: 


^(Ar»-V')' = 0 

lA' 1 a'Y /A' A-1\/1A' 1 _ 

2A c?+l(j/”*^\A~*' r / \2A d + 1 a) 


(9) 

( 10 ) 


[2^ 



+ (iV-l) 



4 


D-2 




/A' D-2 a'\ 

“ {d+l)iN-2)^) 

1 cr'Y fl 1A'\ 
N-2~^ / 2X ) 


flA' a'\ 1 d A'a' 

\2 A~ a ) ^ 2d+l A a 


fl 1 ( t '\ 

\r^ N -2 a ) 


( 11 ) 


A' 


A' 


A-2A-1 


1 


N-2 


N-l\ 

a' 

= 0 

r ) 

a 


(A' 

N 

1 _ 

V A 

n- 


= 0 


( 12 ) 

(13) 


where the prime denotes derivation with respect to r. 

The third group contains one equation, which involves both the time deriva¬ 
tive and the radial derivative of the field variables: 


SfN-lA' D-2a'\ dt A' 4 ■ 

S \ 2 ~A ~ d+l~a ) ^ 2T~A ~ D-2^^ ' 


(14) 


This equation (14) comes from the (rt)-component of the Einstein equations. 
The equations (9-13) have a solution 



(15) 

a{r) = {A{r)y 

(16) 

1/p oA 7 ( 1 - 7 ) A' 

2^ ’y {d+l){N-2) A 

(17) 


where tq is an integration constant, and the constant 7 can take any value 
in the range 0 < 7 < 1, as long as only the equations (9-13) are taken into 
consideration. 

On the other hand, the equations for the scale factors can be solved by 
assuming a power law behaviour of the scale factors such as 


S(xe and T oc 


(18) 
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where a and f3 are constants. Furthermore, if we take the time derivative of the 
scalar field as 

(19) 

where <5 is a constant, then the equations (5-8) yield the following relations 
among a, (3 and 5: 

Na + d^=l (20) 

Na ^ + d/32 ^ ^ j ^ (21) 

Apparently, these relations are a generalization of the Kasner condition in 
higher dimensions. 

Finally, the equation (14) gives a relation among a, /3, 7 and <5. We find 5 
can be solved in terms of a and 7 and obtain 


1 j {d+l){N-2) ( 1-a D 

7 ( 1 - 7 ) I 2 d 



where the sign in the right-hand side matches that in equation (17). The scalar 
held conhguration can be expressed, for instance, as 


.±2(^-2) 


7(1 - 7) 


{d+l){N-2.) 
{d+l){N-2) ( 1-a 


ln(A) 


7(1-7) \2{D-2) 


(23) 


where (/)o and to are integration constants. 

Since there are three equations (20), (21) and (22), there is only one degree 
of freedom in choosing the values of the constants. For a given value for a 
constant, say 7 (0 < 7 < 1 ), the values of the other constants are given by roots 
of quadratic equations in general. 

Now we examine two special cases, for d = 0 and for N = 3. 

Case 1. d = 0 {D = N + 1). In this case the metric does not include the 
extra scale factor (T); therefore the values for the constants a, 7 and 6, which 
appear in general cases, are hxed as: 


1 

N 

(24) 

1 / N 

(25) 

^ “ 2 ^ Y 8(Af- 1) 

N- 1 

(26) 

d = ± ^ 

2VN 


where the sign in equation (26) is independent of the sign in (25). 
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Now the spherically symmetric metric takes the form 




dr"^ 

A(r 


+ r^dn"^ 


N-l 


where 


aw = i-(I5) 


N-2 


1 


a(r) = (A(r))’^ with 7 = — ± 


N 


2 V8iN-l) 


S oc . 


The scalar field is then 

1 I II 


= d’n ± 


{N - 1) ln(A) + sign ( 2 “ 7 


1 \ N-l rt 

In 


yiv V^o 


(27) 

(28) 

(29) 

(30) 

(31) 


Case 2. N = 3 {D = d + 4). In this case the value for the constants a, /3, 7 
and 6 are determined by solving the equation for a and 7 


3a^ 


(1 — 3a)^ 


1 


dT2 / d-t-11 — cx 


7(1 — 7 )c?+l\fi + 2 2 


07 


= 1 


(32) 


and the other equations yield the values for jd and 5. The contours for possible 
values for a, and 7 are plotted in figure 1 , where the contours for d = 1 and 
for d = 00 are shown. As seen from figure 1, we can find a set of parameters 
leaving one degree of freedom for d ^ 0 , in general. 

In the following section we consider the exact solution to the system governed 
by the action (1) with non-zero A. 


3 Solutions for a scalar field with an exponential 
potential 

Now we turn to the A 0 case. The held equations then take die following 
form: 

n<j3 = (33) 

Rmn = m 4N/N(l> + 9mn ■ (34) 

We use the same ansatze for the metric and the dilaton scalar held as in 
the previous section. To obtain analytic solutions we assume that the equations 
(9-14) are unchanged even if there is a potential term. This assumption with 
the homogeneous cosmological solution [ 10 ]. 
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Figure 1: The contours for possible values for a. and 7 in the exact solutions 
obtained for a massless scalar field and = 3 (in section 2); the contours are 
for d = 1 and for d = 00 , as indicated. 


We find that the assumption requires two constraints on the field variables 
as follows: 

The differential equations including the scale factor, which correspond to 
(5-8), are reduced to: 






,4a0/(D-2) 


(36) 






A 


D-2 
A 


cr2/(d+l) 

A 


„4a0/(D-2) 


The equation (15) becomes 

P-2)| 


1 A' 

2 ^ 


D-2 \a^/(d+i) 


1 cr' 


4a<j>l(D-2) 


d + 1 a 


D-2 


(37) 

(38) 


(39) 


The differential equations (9-14) remain vaild in this case. Thus they again 
call for the functions of r as 

A(r) = l-(^)'^”' (40) 


7 




cr(r) = (A(r))'>' 


( 41 ) 


/ = ±i / 7(1 - 7) A' 

2 ]] {d+l){N-2) A 


(42) 


where 7 is left undetemined in this step. 

As long as the metric (4) is used, the solutions for S are treated separately 
for a ^ 0 and for a = 0 . 

For a 0 we have the following solutions to (35-38) 


A oc 


(43) 


D-21 


2a t 


with 


A 


.4a0/(D-2)_ (^-2)(i?-l-a2) 1 

cr2/id+i) Aa4 t2 


(44) 


(45) 


Note that the sign of the time derivative of <j) is given definitely in this case 
(A ^ 0). 

Substituting (40-44) into (39), we find that 7 is a root of the following 
quadratic equation 


a2(d+l)7(l-7) fl 


N-2 

and the scalar field is expressed as 


= ( 


(46) 


D-2 

4a 




where (j)Q and to are constants that are mutually related through the relation 
(45). The equation (46) has real solutions if and only if > {N — 2){d — 1). 
This constrains the possible range of a^ if d > 2. 

Turning to the case with a = 0, we have the solutions to (35-38): 

S' oc (48) 


with 


= 


A 


and 


{D-l){D-2) 
6 = Q. 


In addition it is required that 


7 = 1) 


(49) 

(50) 


(51) 


which comes from (35). For a = 0 we obtain the exact solution only for d = 0 
and for d = 1 . 
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(52) 


• d = 0. Then 7 = 5 - The scalar field is expressed as 

(j) = (j)o± — . ln(A) 

where 0 o is a constant. 

The metric then takes the form: 

. (53) 

J 

• d = 1. Then 7 = 1 . The scalar field is constant everywhere, 

(j) = (j)o. (54) 

The metric then looks like 

ds^ = -dt^+ + I . (55) 

The global structure of the spacetime expressed by the exact solutions will 
be examined in the subsequent section. 


4 The global structure of the spacetime 


We examine the global property of the spacetime described by the solutions 
obtained in sections 2 and 3. Here we take the d-dimensional space as an 
extra space, and so we study the {N -\- l)-dimensional spacetime as a physical 
spacetime. Therefore we will concentrate attention on the time and radial part 
of the metric. Furthermore, it is useful to transform the time coordinate t into 
7], which satisfies 

■ ( 56 ) 


' s^{t) ■ 

Then the ( 77 , r) part of the metric is expressed as 


ds^ = S^ir)) -CT-2/('^+i)(r)A(r)d772 + 


^2/{N 2)(r) 2 


Here we must recall 

A(r) = 1 — ) and a(r) = (A(r))^ (58) 

where tq is an arbitrary constant. 

Using this ‘conformally invariant time coordinate’ 77 , the scale factor S in 
the case with A = 0, treated in section 2 , can be written as 

S{tj) = {Arj + (59) 
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where A and B is arbitrary constants. We remember that there is a rehtion 
between a ard 7 , which can be obtained from ( 20 - 22 ), 


iVa^ + 


(l-fVa)2 

d 


(d+l)(iV- 2 ) 

P-2)7(1-7) 


(1(1-a) 


D-2 Y 


= 1 . 


(60) 


For the A 0 treated in section 3, the relations among the several constants 
take different forms. In the coordinate system that includes the conformally 
invariant time, the expressions for the solutions take dissimilar forms for a 1 
and for a = 1 . 

For a 1, the solutions can be written as 


5 ( 7 ) = + 


(61) 


with 


D-2 

4a 


1 - 


27 

d + 1 


2(1 - 1 ^ 2 ” 


a^(d+ 1)7(1 - j) _ fl 
N-2 


= ( 2^^+ “"I 


(0 < 7 < 1) 


p 4 a 0 o/(£>- 2 ) ^ (£> 2){D 1 |2 


A(l-a 2)2 


(62) 

(63) 

(64) 


where A, B and (po are constants. The value for must be greater than 
{N — 2){d — 1), for the equation (63) has real roots. 

For a = 1 we find the form: 


S{r)) = (65) 

<P = (po- ^{D - 2) ln(A) + ^{D- 2)h{ri - rjo) (66) 

with 

(Ath2(^ = (i(,+ l)_,)’ ,0<,<1) (67) 

j»./(D-2) ^ lA_2l!fc2 (68) 

where h, rjo and po are constants. Incidentally, the equation (67) has real roots 
only for d = 0,1 or d = 2 and N = 3. 

The structure of the spacetime singularity involved in this exact solution is 
as follows: the timelike singulality is located at r = ro, regardless of the presence 
of A and of the value for a.^ However, the property of another ‘cosmological’ 
singularity depends on the value for A and a. 

Let us consider the contracting universe, i.e. the case with the scale factor 
decreasing with time. This corresponds to choosing A < 0 in (59) and A > 0 in 

®This is guaranteed by the fact that [{D — 2)/{{N — 2){d + 1))]7 < 1 . 
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(61), and ft, < 0 in (65). We further assnme that the values of constants B and 
? 7 o are zero. 

For A = 0 the spacelike singularity is located at r/ = 0, which is just the ‘big 
crunch’ singularity. The range of the cosmological time is — oo < ry < 0. The 
spacelike singularity lies at rj = 0. Inward going radial light rays reach either 
the timelike or the spacelike singularity, while any outward going radial light 
rays reach the null singularity. 

On the other hand, for A 0, the global structure of the singularities de¬ 
pends on the value of a (see figure 2). 




A=0 







AseO 




Figure 2: The schematic view of the global property of the spacetime described 
by the exact solutions obtained in sections 2 and 3. 


For a > 1 the global structure of the singularities is the same as the case of 
A = 0. The spacelike singularity lies at ry = 0 (the range of ry is — oo < ?y < 0). 

For 0 = 1 there is a null singularity at ?y = oo (the range of rj is — oo < rj < oo) 
[15]. Inward going light rays hit the timelike singularity at r = tq, while any 
outward going light rays reach the null singlarity. 

For 0 < o < 1 the cosmological singularity is still null. But the range of the 
cosmic time is 0 < ry < oo. 

For 0 = 0 there is no other singularity besides the timelike one at r = rg, 
because the spacetime is asymptotically de Sitter space in this case. 

Next we examine apparent horizons in these spacetimes. It is significant 
to study the property of apparent horizons, especially when the metric is not 
static. The hiding of singularities by apparent horizons may imply the formation 
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of black holes in many cases, for example. 

An apparent horizon is Ihe surface defined by 


g^^R,MR,N = 0 


( 69 ) 


where 

R{ri,r) = . 


In our case the equation (69) leads to 


_ 1 
S dr] r 




Af-2l 




N-2 


D-2 

■ iM-2){d+i)y 


(70) 


(71) 


This defines the apparent horizon in the spacetime. 

The apparent horizon may be spacelike, null, or timelike in different regions 
in general cases. The property can be indicated by the ratio of the slope of the 
apparent horizon and the null ray [ 8 ]. The absolute value of the ratio is given 
in our case as 


dr]AH 

drjNuii 


a 

1 — a 


D-2 N-2[l- (ro/r)^-"][l + d(l - l){ro/r)^-^] 

d+1 d+1 [1 - (1 - 7 )(ro/r)^“ 2]2 


(72) 

for the case with massless scalar. For the case with the exponential scalar 
potential (a ^ 1) equation (72) is still valid if a in the equation be replaced by 
1/a^. For 0 = 1 equation (71) merely determines the location of the timelike 
apparent horizon. Now, let us examine the value of (72) for each exact solution 
obtained in sections 2 and 3. 

First we examine the massless scalar case, treated in section 2. For d = 0, 
i.e. thee is no ‘extra’ dimension other than the spherically symmetric space, the 
equation (72) then becomes 


drjAH _ ^ N-2 1 - (rp/r)^ ^ 

drjNuU iV - 1 1 - (1 - 7)(ro/r)^“2 


(for d = 0, massless) (73) 


since the value of a, has been uniquely solved as a = 1/N (24). In this case 
die apparent horizon is spacelike for r > tq. Husain et al found this feature 
for the N = 3 case [ 8 ], and we find there that the same characteristic of the 
apparent horizon holds for arbitrary dimensions of spherically symmetric space. 
For d > 1 the apparent horizon may be spacelike, null, or timelike in different 
regions, in general. We can, however, and a set or parameters (a and 7 ) which 
allows a spacelike apparent horizon in all regions. One cannot choose a set of 
the pammeters that admits a timelike apparent hoizon in all of the spacetime 
region, except for = 3 and d> 2. 

Next we turn to the case with A ^ 0. The property of the apparent horizon 
depends on a as well as N and d. It is classified into four cases: 


• a = 0. To require an exact solution of the type we have considered in this 
paper, the allowed values for d are d = 0 and d = 1. For each case the 
apparent horizon is timelike, except at i] = 0, r = tq. 
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• 0 < o < 1. In this case d = 0 and d = 1 are permitted as well. The 
apparent horizon may be timelike, null, or spacelike in a different region. 
But in the viciniy of r = rg and at spatial infinity the apparent horizon is 
always timelike. 

• a = 1. The allowed values for d are d = 0, d = 1, and d = 2. Further, 
for d = 2 only = 3 is pemitted. For each case the apparent horizon is 
timelike everywhere in the spacetime. 

• a > 1. In this case any value for d is permitted, but the value for N is 
restricted by a and d for d > 2 through the inequality > (A^ — 2)(d— 1). 
The apparent horizon may be timelike, null, or spacelike in a different 
region. For a sufficiently large a the apparent horizon becomes spacelike 
everywhere. 

The schematic view of the global structue of the spacetime represented by 
our exact solutions is exhibited in figure 2. The cases in which the spacelike 
apparent horizon covers the future singularity can be regarded as models of a 
gravitational collapse, provided that the timelike singularity can be ignored. 

5 Masses for the spherical system 

In this section we evaluate the mass of the self-gravitating system of a scalar 
field described by the exact solutions obtained in section 2. 

In the present analysis we adopt the exact solutions with d = 0, which de¬ 
scribes spherically symmetric spacetimep in this case the exact solution can be 
regarded as a model of gravitational collapse, since the future spacelike singu¬ 
larity is covered by the apparent horizon. The mass can be defined by local 
variables in the spherically symmetric system [16]: 

m{r], r) = R,m R,n ) (74) 

Idtt 

where Am-i = 27r^/^/F(Af/2) and R is defined by (70). 

On the apparent horizon, the mass can be written in the following form: 

~ 167 r [I-(l-7 )(ro/r)^-2](^-2)/(Af-i) • 

(75) 

This expression exhibits a complicated form, which depends non-trivially on the 
spacetime dimension. 

The exact solution obtained in section 2 may not be appropriate for a model 
of gravitational collapse, because of the singular behaviour of the scalar field 
configuration as an initial condition. Our result, however, suggests that some 
physical quantities may depend on the space dimension. Thus we also feet 
interest in numerical study of the higher-dimensional gravitating system.^ 

^The authors of [8, 12, 13, 17, 18] are also interested in the analytical interpretation of the 
numerical results. 
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6 Conclusion 


In this paper we have given exact solutions for a spherical collapse driven by 
a scalar field with an exponential potential in Z?-dimensional spacetime. The 
solutions describe evolution of a scalar field configuration in the background 
metric of a product of a spherically symmetric space and an internal space. This 
solution could be generalized for various supergravity models, which include 
dilaton fields. 
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